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1.  Introduction 

Consider,  for  example,  a  classical  mechanical  system  with  Lagrangian 

(1.1)  L(x,  i)  =  I  -  U(x). 

The  wave  function  of  the  quantum  mechanical  system  corresponding  to  thig 
classical  one  changes  with  time  t  according  to  the  Schrodinger  equation 

(1-2)  «0+,  .)  =  .(.), 

Feynman  [3]  expressed  this  wave  function  x)  in  the  following  integral 
form,  which  we  shall  here  call  the  Feynman  integral 

(1.3)  \l/{t,x)  =  j  [^“  Uixr)'^dTj  <p{x^)YldXr, 

where  F*  is  the  space  of  paths  X  =  (xt,  0  <  t  S  t)  with  xo  =  x,  Hr  dxr  is  a 
uniform  measure  on  and  iV  is  a  normalization  factor.  It  should  be  noted 

that  the  integral  [^r/2  —  U{Xr)]  dr  is  the  classical  action  integral  along  the 

path  X.  (This  idea  goes  back  to  Dirac  [1].)  It  is  easy  to  see  that  (1.3)  solves 

(1.2)  unless  we  require  mathematical  rigor.  It  is  our  purpose  to  define  the  gener¬ 
alized  measure  Ylrdxr/N,  that  is,  the  integral  /  F(X)  Hr  dx^AT,  rigorously  and 

JTz 

to  prove  that  (1.3)  solves  (1.2)  in  case  U{x)  =0  (case  of  no  force)  or  U{x)  =  x 
(case  of  constant  force).  See  theorem  5.2  and  theorem  5.3  below.  We  hope  this 
fact  will  be  proved  for  a  general  U{x)  with  some  appropriate  regularity  condi¬ 
tions. 

Our  definition  is  also  applicable  to  the  Wiener  integral]  namely,  using  it,  w€ 
shall  prove  that  the  solution  of  the  heat  equation 

(1-4)  If  “  III  “ 

is  given  by 

(1.5)  u(t,  ^  j  exp  U{xr)~^dT^f{xt)  H  dxr 
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for  any  bounded  continuous  function  U{x).  See  theorem  4.3.  This  should  be 
called  the  Feynman  version  of  Kac’s  theorem  that 

(1.50  «((,  x)  =  exp  [- U{x,)  dr]  fix,)W,(dX) 

solves  (1.4).  In  this  paper  the  paths,  that  is,  the  points  in  r^,  are  denoted  with 
capital  letters  X,  Y,  ■  •  •  and  their  values  at  time  r  are  denoted  with  the  corre¬ 
sponding  small  letters  with  the  suffix  r  such  as  Xr,  Vr,  •  •  • .  Now  that  Kac’s 
theorem  is  well  known  to  probabilists,  no  one  bothers  with  its  Feynman  version. 
However,  it  is  interesting  that  Kac  had  the  Feynman  version  (1.5)  in  mind  and 
formulated  it  as  (1.5')  to  make  it  rigorous  [5]. 

Gelfand  and  Yaglom  [4]  proposed  a  method  of  defining  the  Feynman  integral. 
They  replaced  h  with  h  —  ia  in  (1.3)  to  reduce  the  Feynman  integral  to  the 
Wiener  integral  and  defined  the  Feynman  integral  as  a  limit  of  the  Wiener  inte¬ 
gral  by  letting  o- 1  0.  Our  method  is  different  from  theirs  in  the  point  that  we 
define  JXr  dXr/N  directly  and  treat  both  the  Feynman  integral  and  the  Wiener 
integral  on  the  same  level. 

2.  The  mathematical  meaning  of  JJt  dx^/N 

What  Feynman  had  in  mind  for  JJr  dxr  must  be  a  uniform  measure  on 
Rigorously  speaking,  this  measure  does  not  exist.  Therefore,  we  should  define  it 
as  an  ideal  limit  of  a  sequence  of  measures  on  In  order  to  be  able  to  com¬ 
pute  the  integral  (1.3)  or  (1.5),  the  approximating  measures  should  be  concen¬ 
trated  on  the  set  L*  of  all  X  =  (xr,  0  <  r  ^  t)  E  satisfying 

(L.l)  Xr  is  absolutely  continuous  in  r, 

(L.2)  Xr  =  dxrldr  G  L^{0,  t], 

(L.3)  limrioa^r  =  x. 

We  shall  now  construct  a  sequence  of  probability  measures  on  L*  whose 
ideal  limit  is  the  uniform  distribution  on  Let  p(t,  a),  with  c,  r  G  (0,  t],  be 
strictly  positive  definite  and  continuous;  for  example,  p(r,  a)  =  exp  (  — |t  —  o-|). 
Let  ^t(w),  w  E  fl(B,  P),  be  a  Gaussian  process  with 

(2.1)  Ei^r)  =  0,  P(^?,)  =  p(r,  <r). 

It  is  well  known  that  such  a  Gaussian  process  exists.  Since  the  continuity  of 
p(t,  <r)  implies  the  continuity  of  in  the  mean,  there  exists  a  measurable  version 
[2]  of  |r.  Denote  that  version  with  the  same  symbol 
Noting  that 

E  ^r  dr)  =  p(t,  t)  dr  <  -f « , 
p{/„‘«5rfT<+»|  =  i. 


(2.2) 

we  can  see  that 
(2.3) 
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Put 


(2.4)  ==  X  +  n  dS, 

Then  x^^\  for  0  <  t  ^  is  also  a  Gaussian  process  with 


(2.5) 


=  a;, 

E{[x^J^'>  -  a:][rc^’*^  -  x]}  =  w®  Jj  Jj  p{di,  dz)  ddi  ddz. 


0  <T 


Denote  by  the  probability  distribution  of  the  sample  function  of  the 
process  Xr"\  with  0  <  t  ^  Then  is  concentrated  on  L*  and  any  finite 
dimensional  marginal  distribution  of  Pn^,  say  over  coordinates  n,  n,  •  •  • ,  Tm,  is 
Gaussian  with  the  density 


(2.6) 


^1/2  r _ 1  "* 


where  the  matrix  (6,7)  is  the  inverse  of  the  matrix  (va)  with 

(2.7)  Vij  =  IJ'  jj*  p{di,  di)  ddi  dBi,  i  «=  1,  2,  •  •  • ,  m 

and  h  is  the  determinant  of  (6,7).  The  existence  of  (6,7),  that  is,  the  nonvanishing 
of  the  determinant  of  (y.y)  results  from  the  assumption  that  p(t,  a)  is  strictly 
positive  definite. 

Since  the  Gaussian  distribution  (2.6)  tends  to  a  uniform  distribution  on  the 
m-space  in  the  sense  that,  for  any  almost  periodic  function  f{xi,  X2,  •  •  •,  Xm), 

(2.8) 

j  f{xi,  •  •  • ,  a;  J  ^2^2)m7'2  ^xp  hj{xi  -  x)  {xj  -  x)  j  da;i  •  •  •  dx„, 

tends  to  the  Bohr  mean  91l(/)  of  /  as  n  ,  it  is  reasonable  to  say  that  P»  \  for 
n  =  1,  2,  •  •  •  approximates  the  uniform  distribution  on  and  that  JJt  dxj  is 
an  ideal  limit  of  this  sequence. 

N must  be  also  an  ideal  limit  of  a  sequence  of  numbers  {Nn}  such  that  Pn^/Nn 
tends  to  I][t  dXr/N  in  some  sense. 

Keeping  these  heuristic  considerations  in  mind,  we  shall  give  a  mathematical 

meaning  to  JJr  dXr/N,  that  is,  to  the  linear  functional  I  (P)  =  I  Fwn-  dx,/N. 

There  are  many  ways  of  defining  this  functional  in  accordance  with  the  choice 
of  the  sequence  {Nn}.  We  shall  express  I{F)  as  7(P,  Nn)  referring  to  the  sequence 

w. 

Definition. 

(2.9)  I(F,  JV.)  =  lira  ^  f  F(X)Pi’\dX). 

n-»«  iV„ 

The  domain  J)(Ar„)  of  this  functional  I(F,  Nn)  is  the  set  of  all  F  for  which  the 
limit  in  (2.9)  exists  and  is  finite. 
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Fixing  {Nn} ,  we  shall  write  3D(iV„)  as  3l)  and  I (F,  Ff„)  as 
(2.10)  fF{X)]ldx,. 

We  shall  mention  three  interesting  cases. 

(i)  Uniform  integral.  iVn  =  1,  with  n  =  1,  2,  •  •  •.  If  F{X)  is  of  the  form 
f(xri,  x-n,  •  • Xr„)  with  an  almost  periodic  function /(.ri,  .To,  •  •  • ,  Xm),  then  F  G  2) 
and 


(2.11) 


where  2ni(/)  is  the  Bohr  mean  of  /. 

(ii)  Wiener  integral.  Nn  =  l/IIi-Cl  +  n  =  I,  2,  ■  ■  ■,  where  will  be 

defined  in  the  next  section.  We  shall  discuss  the  Wiener  integral  in  section  4. 

(iii)  Feynman  integral.  Nn  —  l/II^Cl  +  n~Xjhiy'~,  with  n  —  1,2,  •  •  •,  with 
the  same  X^  as  in  (ii).  This  will  be  discussed  in  section  5. 


3.  Orthogonalization  method 

In  the  following  sections  we  shall  be  faced  with  the  integrals  of  the  type 

(3.1)  /  =  f  G{X)Pl'\dX). 

J  l^x 

Recalling  that  was  defined  as  the  probability  distribution  of  the  sample 
path  of  the  process 

(3.2)  T(”>(a))  =  .f  +  n  Uco)  de 

introduced  in  section  2,  the  integral  I  can  be  expressed  as  the  mean  value  of 
G[X("Kc«^)]  on  S2(B,  P) 

(3.3)  I  =  f^G[X^'‘>M]I’(<kc)  =  £{G[X(«i]}. 

To  compute  this,  we  shall  use  the  usual  orthogonalization  method.  The  idea  is 
as  follows.  Let  T  denote  the  operator  from  JJ{0,  t]  into  itself, 

(3.4)  ('Ft;)(t)  =  p(t,  cr)v(a)  da. 

Then  T  is  a  strictly  positive-definite  compact  operator.  Therefore  T  has  positive 
eigenvalues  {X„}  whose  eigenfunctions  {r/„}  constitute  a  complete  orthonormal 
system  in  t]. 

Now  put 

(3.5)  a,{o})  =  =  Jp*  ^r(w)i7.(r)  dr] 

this  inner  product  is  well  defined,  thanks  to  (2.3).  Then  {a„}  is  a  Gaussian 
system,  since  is  a  Gaussian  process.  Equation  (2.1)  implies  that 
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(3.6)  Eia^a^)  =  p(t,  (T)nyiT)v,ii<r)  dr  da. 

Therefore  a,,  with  =  1,  2,  •  •  •,  are  independent  and  each  a„  is  Gaussian  with 
mean  0  and  variance  X,.  Since  we  have 

(3.7)  2  X.  =  £  E(ai)  =  £(E  al)  =  E  (^f‘  {f  dr)  =  p(r,  r)  dr, 

the  continuity  of  p(t,  a)  implies 

(3.8)  L  X.  <  +«  ; 

this  fact  will  be  useful  in  the  following  sections. 

Noting  that 

(3.9)  |T(to)  =  a,((a)'ny{T) 

P 

and  that 

(3.10)  4'‘^(w)  =  X  -\-n  ^e{o))  dd. 

we  can  express  I  in  the  form 

(3.11)  /  =  E{H{ax,a^,  •••)} 

with  some  H.  Using  the  independence  and  the  normality  of  {a^} ,  we  can  compute 

(3.11)  more  easily  than  the  original  form  (3.1). 


4.  Wiener  integral 

We  shall  now  discuss  the  integral  (2.10)  for 


(4.1) 


Nn  = 


1 

n  (1  + 

¥ 


n  =  1,2, 


where  X^,  with  v  =  1,2,  •  •  •,  are  the  eigenvalues  introduced  in  section  3. 

We  can  verify  easily  the  convergence  of  the  infinite  sums  and  infinite  products 
appearing  in  this  section  by  appealing  to  (3.8). 

Lemma  4.1. 

(4.2)  Q'«{dX)  ^Y,^xp(^-0dr')pi‘>(dX) 

is  a  probability  distribution  on  Lx. 

Proof.  Using  the  orthogonalization  method,  we  get 

(4.3)  h  =  exp  (  -  I  dr)  P“>  (dX) 

=  p[exp(-«^0dr)] 

=  £[exp(-,i*E|)]. 
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Noting  that  the  a„  ioT  v  —  1,  2,  •  •  •,  are  independent,  we  have 


(4.4) 


/.  =  n£[exp(-^)] 

=  n — — 

V  (1  +  n=X,)'« 


n^a^\ 

2  ) 


da 


=  Nn, 

which  proves  our  lemma. 

Lemma  4.2.  For  any  g  G  /],  we  have 

where  g,  =  (g,  rj,)  =  giT)v.(T)  dr,  for  v  =  1,2,  •  •  • . 

Proof.  By  the  same  idea  as  in  lemma  1,  we  have 

(4.6)  In  =  exp  Xrgir)  drj  Q'^\dX) 

=  i ^  exp [i  J‘ x^(r)  rf’-  -  r I *] 

=  ^j[exp(inZ9^a,-n-zf)] 

=  jl^n^^l^exp  -  n?  |^)  j 

_ 1_TT  f  1  /  ~q:^  ,  .  w^\  , 

(27rX,)»/2  (  2X,  2 

_  1  TT  1  r  ~] 

N„¥(1  +  n%)‘/2  l2{n^X,  +  1)  J’ 

which  proves  (4.5)  by  virtue  of  (4.1). 

As  an  immediate  result  from  lemma  4.2,  we  obtain  the  following  lemma, 
noting  that  J^^gl  =  girY  dr  and  that  gir)  dBir)  is  Gaussian  distributed 

with  mean  0  and  variance  girY  dr  for  the  Brownian  motion  B(t). 

Lemma  4.3.  For  any  g  G  L’'(0,  t\, 


and 

(4.8) 
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^  J"  Gxp  J  XtQ{t)  dr  J  2  dr  J  JX  dXj 

^  ju  Jo 

where  TF*  is  the  probability  measure  for  the  Brownian  motion  process  starting 
at  X,  namely  the  Wiener  measure. 

Theorem  4.1. 

(4.9)  i  exp  (  -  jT'  I  dr)  n  =  WMX) ; 

rigorously  speaking,  we  have,  for  any  continuous  bounded  tame  function  F{X) , 

(4. 10)  F(X)  exp(^-  j'^dr^e  a 

(4.11)  F(X)exp(-  r|dT)dn^,  =  [  F(X)W.(dX). 

A  tame  function  is  a  function  defined  on  an  infinite  dimensional  space  which 
depends  only  on  a  finite  number  of  coordinates. 

Proof.  F(X)  can  be  expressed  as 

(4.12)  F{X)  =  fiXn,  Xn,  •  •  •,  Xt^),  0  <  T1  <  •  •  *  <  Tm  S  t, 

with  a  continuous  bounded  function  /  of  m  real  variables.  To  obtain  theorem  4.1, 
it  is  enough  to  prove  that 

(4.13)  lim  I  f{xru  Xr,,  • .  • ,  XrMn\dX)  =  i  /(Xn,  ,  XrJW.idX). 

n^oo  J^*  yr» 

Let  Qn^  and  Wx  denote  the  marginal  distributions  of  Qn^  and  TF*  over  the  co¬ 
ordinates  Ti,  t2,  •  •  •,  Tn  respectively.  Then 

(4.14)  =  j  ...  j  exp  [i{ziai  -f  •  •  •  -|-  Zma„)]Qij^\dai  •  •  •  dam) 

=  exp  [i{ZiXr,  +  •  •  •  +  Z„XTj)]Qn\dX) 

=  exp  [i{zi  -f  •  •  •  +  2m)x]  j^exp  j  ^(r)xrdrj  Q^n\dX), 

where  gir)  =  5Zr-i2;Y’j(r)  and  ^,(t)  is  the  indicator  function  of  the  set  (0,  tj). 
Using  lemma  4.3,  we  get 

(4.15)  In  exp  [i{z^  +  ...  +zm)x  +  i  I  g(r)  dx.]TF,(dX) 

=  /r,  [iiZxXr,  -!-•••+  ZmXr^]Wx{dX) 

~  I  I  [^(^i“i  +  '  •  '  +  ZmOtrf)]W x{dai  da2  •  •  •  dan). 

Therefore  Qn^  — >  Wx  in  the  weak  sense  as  n  — » <» ,  which  implies  (4.13). 
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Theorem  4.2. 

If  f{x):  C  is  continuous  and  hounded,  then 

(4.16) 

exp(-^  ^dr'^fixt)  E 

and 

(4.17) 

u{t,  =  ^  j 

solves 

(4.18) 

|  =  .(0+,.)=/(x). 

Proof.  Using  the  previous  theorem,  we  have 

(4.19)  «((.  x)  =  f  f(x,)W.(dX)  =  j  fiv)  exp  dy 

and  this  solves  (4.18). 

Theorem  4.3.  {Feynman’ s  version  of  Kac’s  theorem.)  If  f{x)  and  U (x)  are 
continuous  and  hounded,  then 

(4.20)  exp|-y  +  U{xr)~^dT^f{xt)  E  SD 

and 

(4.21)  u{t,  ^  j  exp|- J  1^1^  +  t/(a!,)  jdr|/(a:0  Jldxr 
solves 


(4.22)  ^  2  aJ  ”  U(x)u,  w(0+,  x)  =  fix). 
Proof.  It  is  enough,  by  virtue  of  Kac’s  theorem,  to  prove  that 

(4.23)  j  ^  ^  U{xr)~^dT^f{xt)Pn{dX) 


j  exp  —j  U{xr)dT  f{xt)WxidX). 


Denoting  the  integrals  in  (4.23)  by  /„  and  /  respectively,  we  obtain 


(4.24) 

In  =  exp  U{xr)  dTjf{xt)QU\dX) 

where 

(4.25) 

I nm  "f“  Rnmy 

=  i;  f  ■■■  [‘  [  U{Xr)  ■  ■  ■  U(x.,)Qi.‘'(dX)  A,  • 

*'=1  *'•  Jo  Jo  JLx 

d/T  vj 

r4.26) 

where  |1 
n  as  m  — 

lU 
>  00 . 

means  the  uniform  norm.  Therefore  tends  to  /„  uniformly  in 

Using  theorem  4.1,  we  have 
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(4.27) 


and  it  is  easy  to  see  that  /«,„  — »  /  as  n  — >  «> .  Taking  the  uniform  convergence 
of  limTO-400  Inm  =  In  ioto  account,  we  have 

(4.28)  lim  In  =  lim  lim  Inm  =  Hm  lim  Inm  =  lim  /«,m  =  I, 


n-*»oo  m— >00 


f»— >00  n— >00 


which  completes  our  proof. 

5.  Feynman  integral 

In  this  section  we  shall  discuss  (2.10)  for 

1 


(5.1) 


Nn  = 


n(-ttr 


n  =  1,  2, 


As  in  section  4,  we  can  easily  verify  the  convergence  of  the  infinite  sums  and 
infinite  products,  by  appealing  to  (3.8). 

Lemma  5.1.  If  Re(6)  >  0  and  c  is  real,  then 

(5.2)  j  exp  (—ha^  +  ic)  da  —  (””4b)’ 

Proof.  This  is  true  if  6  >  0.  By  analytic  continuation,  we  can  verify  (5.2) 
for  Re(6)  >  0. 

Lemma  5.2.  If  gr(r)  G  L^iO,  t],  then 

(5.3)  i  f  exp  [I  r  I  dr  +  i  r  x^(t)  dr]  n”>(dX) 

n  «  n^\vhtQy  ~\ 


where 

(5.4) 


Q*  =  {Q,  Vv)  =  ^  g{r)Vy{r)  dr. 


Proof.  We  shall  use  the  orthogonalization  method  introduced  in  section  3. 

"t 


(5.5) 


=  E  [exp  (^  al  +  inp,o.)] 


2rv2 


+ 


in‘'a 

2h 


+  ing, 


da. 
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Using  lemma  5. 1  to  evaluate  this  integral,  we  have 


7„  =  exp  -  X; 


2(n%  -{-hi). 


Noting  that  girY  dr,  we  obtain,  as  an  immediate  result  from  lemma 


Theorem  5.1.  If  gir)  G  IY{0,  t\,  then 


(5.7)  exp  j  ^  dr  +  z  j  Xrgir)  dr  j  G  © 
and 

(5.8)  ^  j  expj^l  j  ^dr  +  i  j  iTg(r)  dr  =  ^xp 

Theorem  5.2.  If  the  Fourier  transform  of  <p{x)  is  a  continuous  function  with 
compact  support,  then 

(5.9)  exp  f  ^dr)  <p(Xr)  G  3^, 


(5.10) 


(5.11) 


exp  Q  J  ^ 

I'it,  ^  f  ^^^(h  fo^  dXr 


i  dt  2  dx^’ 


i/'(0,  x)  =  (p{x). 


Proof.  It  is  enough  to  prove  that 


(5.12)  r  ^  *) 


tends  to 


—  (nr.  — 


<p{y)  dy. 


/  {2Thit)^i^  _  2hit  _ 

Denoting  the  Fourier  transform  of  (p{x)  by  (p{(t)  or  (5F(p)  {x)  as 
(5.14)  <pi£)  =  (?F^)(f)  =  /:.  exp  {2‘7nJtx)(p{x)  dx, 


(5.14) 


we  have 


(5.15)  In  = 


=  Y  I  /*  exp(-27rtfx,).^(a!)diP^^HdX) 

= -^  J  exp  (—2'!ri(tx)  dA  j  exp  j  ^  dr 

—  2TdA  j  XrdT^Pn^d 
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Putting  g{i)  =  —2TrxXr  in  lemma  5.2  to  compute  this  integral  over  Lx,  we  have 


/„  =  y*  (p{A)  exp  j^— 27rtia;  —  ^ 


-:r  dx. 
,2)J 


(5.16)  . „  -  ,  ^  2(^2X.  +  hi). 

Recalling  the  assumption  that  ^(i)  has  a  compact  support,  and  noting  that 

(5.17)  £  g‘‘(r)  dr  = 

we  have 

(5.18)  lim  In  =  f  ^(x)  exp  (—2Tnix  —  2Tr^hit0)  dx. 

n — >00  J 

Since  the  Fourier  transform  of 


(5.19) 


N{x,  hit)  = 


1 


75  exp 


\  2hit ) 


{2irhityi^ 

in  the  Schwartz  distribution  sense  [6]  is  exp  {—2ir^hiW),  we  obtain 

(5.20)  lim  In  =  ^^0]} 

n^oo 

=  <p{x)  *  N{x,  hit), 

namely 

(5.21)  \f/{t,  x)  =  N{x,  hit)  *  (p{x), 
which  completes  our  proof. 

Theorem  5.3.  If  the  Fourier  transform  of  (p{x)  has  compact  support,  then 


we  have 

(5.22) 

exp  J  ~  ^  ^ 

and 

(5.23) 

4'{t,x)  =  ^  J  exp  {^-x^dr  ip{xt)WdXr 

solves 

(5.24) 

hdri/  h^  d^I/  ,  /mi  ^  (  \ 

WO+.o:)  =^(x). 

Proof. 

Defining  ^(:6)  as  in  (5.14),  we  obtain 

(5.25)  In 

r  exp[|  r(|-x.)dr]v(x,)P‘,'>((iX) 

f  t^(£)  di  f  exp  (j  r  f  *  -  1  r  as.  dr  -  2^x,)  P?\dX) 

[  m  exp  [  -2«  (i  +  ^)  x]  dd  ^  exp  [|  J‘  |  dr 

+  ij\(r)x,dT  P^\dX), 
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where  gij)  =  —(t  —  T)/h  —  2vx.  Using  lemma  5.2  to  evaluate  this  integral  over 
Lx  and  recalling  that  ^  has  compact  support  to  take  the  limit  as  n  — >  oo ,  we 
have 

(5.26)  x) 

=  I  m  exp  [-2i  (^  +  ^)^  -  I 

Thus  we  get 

(5.27)  Ul,x)  =^[<1,(1,  ■)] 

=  “  2^)  “P  [( +  O'  ”  (-2-i)>]}- 

Simple  computation  shows  that  ^{t,  x)  satisfies 

iff  =2 m’  ^(0+,i)  =  ^(i)- 

This  implies  (5.13). 

I  would  like  to  thank  Professor  G.  Baxter  for  his  valuable  suggestions  and 
assistance  in  writing  the  final  version  of  this  paper. 
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